Charged-particle rapidity density in Au+Au collisions in a quark combination model 
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Rapidity /seudorapidity densities for charged particles and their centrality, rapidity and energy 
dependence in Au+Au collisions at RHIC are studied in a quark combination model. Using a 
Gaussian-type rapidity distribution for constituent quarks as a result of Landau hydrodynamic 
evolution, the data at ^/snn = 130, 200 GeV at various centralities in full pseudorapidity range 
are well described, and the charged particle multiplicity are reproduced as functions of the number 
of participants. The energy dependence of the shape of the dN c h/dn distribution is also described 
at various collision energies y^ijvjv = 200, 130, 62.4 GeV in central collisions with same value of 
parameters except 19.6 GeV. The calculated rapidity distributions and yields for the charged pions 
and kaons in central Au+Au collisions at tJsnn = 200 GeV are compared with experimental data 
of the BRAHMS Collaboration. 
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O ■ I. INTRODUCTION 

£: 

• The relativistic heavy ion collider (RHIC) at Brookhaven National Lab was built to search for quark matter, or the 
so-called quark-gluon plasma (QGP). Since its first run in 2000, a huge number of data have been accumulated and a 
t-H 1 comprehensive analysis of these data has been carried out. A variety of experimental facts from different aspects imply 
that the strongly coupled QGP has probably been produced in central Au+Au collisions at RHIC. For recent reviews 
of QGP and summary of experimental data, see e.g. [l], 0, 0, 0, 0, 0|- Central Au+Au collisions are characterized 
by the production of thousands of charged-particles in vacuum. The charged-particle density per unit rapidity or 
pseudorapidity dN c h/dy or dN c h/drj is one of the most important observables to measure for the signal of QGP, from 
which a lot of information about the hot and dense matter can be extracted @, 0, 0, OH HH, G2> Ell • One can scale 
dN c h/dr] or dN^/dy by the number of participant nucleon pairs (A^ ar t /2) and observe its logarithmic increase with 
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(iVp ar t), which is regarded as an evidence of color glass condensate [8|, [9|, [14|, [15|]. From the rapidity /pseudorapidity 
density and the transverse energy per particle, one can determine via Bjorken method the real density of the fireball 
which can provide one piece of evidence for the deconfinement phase transition. The ex per imental data about 
__!_ ' the charged-particle rapidity density have been presented by the PHOBOS collaboration [lfil [l7| . the PHENIX 
collaboration [Hj], and the BRAHMS Collaboration [H,[20j]. 

In this paper we will use a quark combination model to study the rapidity/pseudorapidity density varied with the 
number of participants and the energy in full rapidity range. The quark combination picture is successful in describing 
many features of multi-particle production in hadronic collisions. In ultra-relativistic heavy ion collisions at RHIC 
energies, a lot of new features are found, e.g. the high ratio of rip/n^ ~ 1 at intermediate transverse momenta, which 
supports quark coalescence or recombination picture [2lL I22I |23| . The quark number scaling of the elliptic flow is also 
a manifestation of the quark coalescence or recombination |24l . l25l . |2(| . In this paper we will use a binary potential 
model for the constituent quark production and then let the constituent quarks combine into initial hadrons according 
to a quark combination rule. Then we allow the resonances in the initial hadrons to further decay to final hadrons 
with the help of the event generator PYTHIA 6.3 [27| . 

The paper is organized as follows. In the next section we give a brief description of the model for constituent quark 
production and combination. In section III, we present our predictions for the rapidity/pseudorapidity densities varied 
with the number of participants in the full rapidity range at ^/JWn — 130, 200 GeV, the energy dependence of the 
dNch/dr] distribution at various collision energies for central collisions, and the results for the rapidity densities dN/dy 
and yields for charged pions and kaons in the central collisions at ^s^n — 200 GeV. The summary and discussions 
are in section IV. 
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II. THE QUARK PRODUCTION AND COMBINATION MODEL 



In this section we give a brief introduction of the quark production and combination model we use. The model was 
first proposed for hi gh e nergy e + e~ and pp collisions [H, [2i| [H HH, [H, HH and recently extended to ultra-relativistic 
heavy ion collisions [3J, |35|. It has also been applied to the multi-parton systems in high energy e + e" annihilations 

USES!!. 
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A. An effective model for quark production 



The quark production from vacuum is a very sophisticated non-perturbative process. The color glass condensate 
model is a semi-classical QCD effective theory for the quark production in heavy ion collisions [bH . [H| . In this paper 
we use a simple model for quark production which is of statistical nature without dynamic details. We determine the 
number of constituent quarks by the total effective energy for producing quarks from the vacuum excitation. The 
effective energy consists of the part for quark static masses and that for effective interquark potentials. 

Consider a system of N q quarks and anti-quarks excited in vacuum , the number of light and strange quarks/anti- 
quarks follow the ratio N u : Nd : N s = 1 : 1 : A s with N q = N u + Nd + N a , where X s < 1 is the strangeness 
suppression factor due to the heavier mass of strange quarks/anti-quarks. The average quark mass is given by 
m = (2m u + A s m s )/(2 + A s ), where m u — md is the light quark mass and m s the strange quark mass. 

We assume that the interaction is characterized by an inter-quark potential V which takes a substantial fraction of 
total effective energy. The constituent quark number can be determined from the following energy equation, 

E = {N q )m + ^-({N q )-l){V), (1) 

which gives the number of constituent quarks as 

(N q )=2{(a 2 +{3E)V 2 ~a], (2) 

where 

Note that the quark number N q follows a specific distribution, so does the potential, we have taken their averages 
in the above equations. In Eq. (|T]), we only included the two-body potential leading to a E x l 2 ~ s 1 / 4 asymptotic 
behavior for N q at high energy if (V) is constant. For a strong coupling system, it is possible that the n-body (n > 2) 
potential might be more important, and the asymptotic behavior then becomes N q ~ s 1 / 2 ™. When n is large, N q 
more and more approaches a logarithmic increase with energy. 

The PHOBOS experiments have shown that above SPS energies, the total multiplicity per participant pair 
(Nch) /{Npaxt/2) in central events scales with ^/snn in the same way as e + e~ collisions [13, [4Cj]. This suggests a 
universal mechanism of particle production mainly controlled by the amount of effective energy available for particle 
production. Based on this property we extend the above quark production model originally applied to e + e~ annihila- 
tion and pp collisions to heavy ion collisions. The average quark number in nucleus- nucleus collisions can be written 
as 

(N q ) = 2[(a 2 + /3^J^) 1/2 - a](N paxt /2), (4) 

where we have taken E = ^/snn- Note that the effective energy E is equal to collision energy y/s for light quark 
events in e + e~ annihilation ( but E ^ \/s for heavy quark events). Here we have assumed that the colliding nuclei 
is fully stopped and all the colliding energy ^snn is used for particle production, which is reasonable for the central 
collisions due to the fully rescatterings of partons towards the local equilibrium. It is also consistent with the Landau's 
space time evolution picture |4lj |. Note that the average number of quarks and antiquarks (N q ) includes not only new 
produced quarks and antiquarks, but also some additional quarks left by the incident nuclei. 



B. Model for quark combination 



In this subsection we briefly summarize how quarks combine into hadrons in our model. Different from many 
coalescence models which do not distinguish directly produced hadrons from final state hadrons, our quark combination 
model only describe the hadronization of initially produced hadrons including resonances. Then all resonances are 
allowed to decay into final state hadrons. Here we make use of the event generator PYTHIA 6.3 [27j to deal with 
resonance decays. The basic idea is to put N q quarks and anti-quarks line up in a one-dimensional order in phase 
space, e.g. in rapidity, and let them combine into initial hadrons one by one following a combination rule. See section 
II of Ref. [HI for short description of such a rule. We note that it is very straightforward to define the combination 
in one dimensional phase space, but it is highly complicated to do it in two or three dimensional phase space [42| . 
The flavor SU(3) symmetry with strangeness suppression in the yields of initially produced hadrons is fulfilled in the 
model [28l. [30j. Using the model, we have described most of multiplicity data for hadrons in electron-positron and 
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TABLE I: Yields of charged pions and kaons compared with BRAHMS data at 200 GeV H||. 
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7T + 


1660 ± 15 ± 133 


1676 


7T~ 


1683 ± 16 ± 135 


1680 


K+ 


286 ± 5 ± 23 


280 


K~ 


242 ± 4 ± 19 


242 




FIG. 1: The centrality-dependent multiplicities of charged particles in Au+Au collisions at ^/snn = 130, 200 GeV. The solid 
lines are our results. The data are taken from PHOBOS [17| . 

proton-proton/anti-proton collisions [28|, [2^, [3(J HH, HH, HH . Also we solved a difficulty facing other quark combination 
models in describing the TASSO data for baryon-antibaryon correlation in electron-positron collisio ns l&Hl- Combined 
with the color flow picture [37[ , the model can describe the hadroniztion of multiparton states [2a, l3 8| . |39l| . We have 
extended the model to reproduce the recent RHIC data for hadron multiplicity ratios, pr spectra [34[ and elliptic 
flows [35j in central rapidity region. 

III. RAPIDITY AND PSEUDORAPIDITY DENSITIES 

In this section, we use our combination model to compute the centrality dependence of distributions of rapid- 
ity/pseudorapidity densities in Au+Au collisions at y/s^N — 130, 200 GeV and study the energy dependence of the 
shape of the dN c h/dr] distribution at various collision energies y / s]vjv = 19.6, 62.4, 130 and 200 GeV for central colli- 
sions, and also calculate the rapidity densities dN/dy and yields for charged pions and kaons in the central collisions 
at ^/IWn = 200 GeV. 

First we have to fix the parameters of the model. There are two parameters m and (V) or a and j3 in Eq.pJ. As 
we pointed out that these quarks and antiquarks are constituent ones, so we use constituent masses m u — = 0.34 
GeV and m s = 0.5 GeV giving the average mass m = 0.36 GeV. The strangeness suppression factor is chosen to 
be A s = 0.55 by fitting the data at RHIC energies [Hj]. The parameter (3 is set to 3.6 GeV -1 which described the 
e + e~ data. The parameters controlling the total multiplicity are the number of quarks and that of anti-quarks. In 
electron-positron and proton-antiproton collisions, the number of quarks is equal to that of anti-quarks, i.e. there are 
no excess baryons in contrast to anti-baryons. For nucleus-nucleus collisions, however, there are some excess baryons 
deposited by the colliding nuclei. The total number of quarks and anti-quarks (N„) is given by Eq.(f4]). The number 
of net quarks can be further determined by the ratio of anti-proton to proton [241 . At 130 and 200 GeV, we find that 
the net quark numbers are about 420 and 360, respectively. 

With these parameters, we calculate the centrality-dependent multiplicities of charged particles at v / s]v7V = 130, 200 
GeV. The results are shown in Fig. [T]and agree with data very well. 

In order to compute the distribution of rapidity/pseudorapidity densities with both energy and centrality depen- 
dence, we have to know the rapidity distribution of quarks and antiquarks before hadronization. 

In the initial state, colliding nuclei are highly Lorentz contracted along the beam direction. After the initial 
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compression phase, the evolution of highly excited, and possibly deconfined, strongly interacting quark matter can be 
described by the ideal relativistic hydrodynamics. Under the assumption of full stopping and isentropic expansion, 
the amount of entropy (dS) contained within the (fluid) rapidity element dy in the Landau hydrodynamic picture is 
given by [H |H MM HI , 



dS 

— = — ttR 2 Iso(3c s exp(/3wf ) 
dy 



(•5) 



where 2/3 = (1 — c 2 )/c 2 and q = ^Juj 2 — c 2 y 2 with c 2 = (^)iscntropic the sound velocity square, cuf is related to 

the initial and freeze-out temperature 7/ and To by u>f = ln(Tf /Tq). R is the radius of the nuclei. 21 is the initial 
longitudinal length, sq is the initial entropy density. Iq and 1\ are the Bessel functions. The quantity itR 2 Isq is fixed 
to normalize the experimental data at mid rapidity. For \wf\ c s y the quantity dS/dy can be approximated by a 
Gaussian distribution, 



dS ^ exp(-^r) 
dy y/2^ 2 



IT (7 



where 



where m p is proton mass, is pion mass, taken Tf rs m„. There is only one parameter c 2 left to be determined 
by experiments. dS/dy is proportional to dN /dy [13] . Therefore, the rapidity distribution for quarks and anti-quark 
before hadronization can be written as 

_ „ cxp(- — ), (8) 

We take the sound velocity square c 2 = 1/4 for QGP before hadronization (c 2 = 1/3 for ideal gas). We have used 
the fact that all quarks and anti-quarks are within the rapidity range y € [—4.2,4.2] at 130 GeV and 200 GeV at all 
centralities of collisions. 

With this input, we give dN c h/dr] as functions of 77 at all available centralities at ^/snn — 130, 200 GeV. The results 
are shown in Fig. [2] and Fig. [3] One can see a good agreement between our model predictions and data in central 
collisions. For peripheral collisions, there is a slight deviation from data. The tails at large pseudorapidities especially 
in peripheral collisions are associated with remnants of collision spectators from the incoming nuclei. Therefore our 
results are slightly lower than data. Now we study the energy dependence of the shape of the 77 distribution of charged 
particles at various energies of Au +Au collisions. We compute pseudorapidity densities in full pseudorapidity range 
in most central collisions at ^snn = 62.4 GeV and compare with BRAHMS data [4jj|. Here, the quark rapidity range 
is also y E [—4.2,4.2] and the sound velocity square is also c 2 — 1/4 same as in 130 and 200 GeV. By studying we 
find that the shape of r\ distribution is mainly determined by the energy and the sound velocity and the the quark 
rapidity region only influences the forward pseudorapidity. The quark rapidity range is possible dependent of collision 
energy but not sensitive to that. The calculation results show that there is same value of the sound velocity for QGP 
in different collision energies above. This indicates a certain kind of universality for the quark matter produced in 
heavy-ion collisions at the late stage of evolution (before hadronization) at collision energies from 62.4 to 200 GeV. 
We also study the results at 19.6 GeV. We find that the model predictions disagree with the data using the constant 
value of the sound velocity square c 2 — 1/4 no matter how the quark rapidity region is chosen. We have to set the 
sound velocity c 2 = 1/7 and quark rapidity range y E [—3.2, 3.2]. The change in the sound velocity might reflect the 
different properties of the matter produced at 19.6 GeV and at 62.4 GeV or higher energies. The model predictions 
are shown in Fig. 2] at ^/snn = 19.6, 62.4, 130, 200 GeV. The agreement with data is also satisfactory, which means 
that our model captures the energy behavior in the available collision energies. 

To separate the trivial kinematic broadening of the dN c h/dr] distribution from the more interesting dynamics, we 
also study the scaled, shifted pseudorapidity distribution dN c h/drf / (N part /2} , where rj' = r\ — yteam , for Au+Au 
collisions at different energies. The calculation results arc shown in Fig. [5] at ^/snn — 19.6, 130, 200 GeV and two 
centrality bins 0-6% and 35-40%. In most central collisions our model can describe the limit fragmentation very well. 
But in peripheral collisions there is a disagreement close to beam rapidity. The reason is that in peripheral collisions 
the beam rapidity region is dominated by spectators and is not covered by our model. 

In ultra-relativistic heavy ion collisions at RHIC energies, charged pions and kaons are copiously produced. The 
yields of these light mesons carry the information on the entropy and strangeness created in the reactions. Here, 
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FIG. 2: Distributions of pseudorapidity density for charged particles in Au+Au collisions at v / sjvjv = 130 GeV for six centrality 
bins. The solid lines are our results. The data are taken from PHOBOS [l7j . 



We calculate the rapidity density dN/dy and yields of charged pions and kaons in full rapidity for central Au+Au 
collisions (0 — 5%) at ^/snn — 200 GeV. The yields of charged pions and kaons compared with BRAHMS data [48[ 
are shown in Tab. I. The results for rapidity density distributions of charged pions and kaons are shown in Fig. [5] 
Here, the pion yields are collected excluding the contributions of hyperon (A) and kaon Kq s decays. One can see that 
our model can well describe rapidity densities dN/dy and yields of charged pions and kaons in the whole rapidity 
range for central Au+Au collisions at ^snn = 200 GeV. 



IV. SUMMARY AND DISCUSSIONS 



We study in a combination model the rapidity and pseudorapidity densities at various collision energies and cen- 
tralizes. We use the Landau relativistic hydrodynamic model to describe the the evolution of highly excited and 
possibly deconfined quark matter. As a result, we obtain the Gaussian-type rapidity spectra of constituent quarks 
before hadronization. Then we use our combination model to describe the hadronization of initially produced hadrons 
including resonances, whose decays are dealt with by the event generator PYTHIA 6.3 [27]. We compute charged 
multiplicities and pseudorapidity densities at a variety of centralities at 130 and 200 GeV. The results for pseudo- 
rapidity densities are in good agreement with data in central collisions. In peripheral collisions, our predictions are 
slightly lower than data due to the fact that our model does not include the influence of the spectators. Our model can 
well describe the dependence of pseudorapidity densities and charged multiplicities on centralities and the number of 
participants respectively. We also calculate pseudorapidity densities at 19.6 and 62.4 GeV which describes the RHIC 
data very well. This means that our model can reproduce the collision energy dependence of pseudorapidity densities. 
However, We find that the value of the sound velocity square c 2 s — 1/7 at 19.6 GeV is different from that at 62.4 
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FIG. 3: Same as Fig. [2] but at ^sJTn = 200 GeV. 




FIG. 4: Pseudorapidity densities dN/d-q for charged particles in most central collisions at various collision energies ^/sjvjv = 
19.6, 62.4, 130, 200 GeV. The lines are our results. The PHOBOS data are from Ref. [H, while the BRAHMS data are from 
Ref. H| 
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FIG. 5: The scaled, shifted pseudo-rapidity rapidity density at y^ijvjv = 19.6, 130 and 200 GeV. The results at two centrality 
bins are presented: 0-6% and 35-40%. The lines are our results. The data are from Ref. [13], 
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FIG. 6: dN/dy in the — 5% most central collisions at y / sjvjv = 200 GeV. The solid lines are 7r" 
are 7r~ and K~ . The data are given by BRAHMS collaborations [48l |. 



and K + , and dashed lines 



GeV or higher energies. This imply that there is large change in properties of the hot and dense matter produced at 
collision energy between 19.6 and 62.4 GeV. To separate the trivial kinematic broadening of the distributions of the 
pseudorapidity density from more interesting dynamics, we compute the scaled and shifted pseudorapidity density 
distributions dN c h/drf / (N par t/2) with rf — r\ — ybeam at collision energies 19.6, 130 and 200 GeV. The good agreement 
with data is found except in the beam rapidity range of peripheral collisions, where our predictions are lower than 
data. Finally we present our results for rapidity densities of charged pions and kaons in most central collisions at 200 
GeV. No contradiction to data is found. Note that the BRAHMS pion data do not include the decay products of K® 
and A, we also make the same corrections. 
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